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ABSTRACT 

B2 1023+25 is an extremely radio-loud quasar at z = 5.3 that was first identified as a likely high-redshift blazar 
candidate in the SDSS+FIRST quasar catalog. Here, we use the Nuclear Spectroscopic Telescope Array ( NuSTAR ) 
to investigate its non-thermal jet emission, whose high-energy component we detected in the hard X-ray energy 
band. The X-ray flux is ~5.5 x I (F 14 erg cm 2 s~' (5-10 keV) and the photon spectral index is T x — 1.3-1. 6. 
Modeling the full spectral energy distribution, we find that the jet is oriented close to the line of sight, with a 
viewing angle of ~3°, and has significant Doppler boosting, with a large bulk Lorentz factor ~ 13, which confirms 
the identification of B2 1023+25 as a blazar. B2 1023+25 is the first object at redshift larger than 5 detected by 
NuSTAR, demonstrating the ability of NuSTAR to investigate the early X-ray universe and to study extremely active 
supermassive black holes located at very high redshift. 

Key words: galaxies: active - quasars: general - quasars: individual (B2 1023 + 25) - X-rays: general 
Online-only material: color figures 


1. INTRODUCTION 

Blazars are radio-loud active galactic nuclei (AGNs) with a 
relativistic jet directed at a small angle from our line of sight 

(Urry & Padovani 1 995). The peculiar orientation relativistically 
boosts the radiation emitted from their jets, making them visible 

even at high redshift (z). 

The typical spectral energy distribution (SED) of a blazar is 

dominated by its non-thermal emission, characterized by two 
broad humps: the lower frequency component is attributed to 

synchrotron emission, while the higher frequency component 
is attributed to inverse Compton (IC) emission. The humps 

of a blazar SED peak at lower frequencies as the bolometric 

luminosity increases, at least for blazars studied to date (see, e.g., 
Fossati et al. 1998; Ghisellini et al. 2011 and Giommi et al. 2012 
for different interpretations of the effect). At z > 4, we expect to 

see only the most powerful objects. Therefore, we should detect 
high-redshift sources whose synchrotron hump peaks in the sub- 

millimeter and whose IC hump peaks in the ~MeV band. The 

synchrotron component shifts far enough out of the optical/UV 

to leave the accretion disk, usually swamped by non-thermal 

emission, visible (Sbarrato et al. 2013; Wu et al. 2013). The 

optical emission of a high-power, high-redshift blazar therefore 
becomes indistinguishable from that of a high-redshift, radio- 

quiet quasar. Since the accretion disk is visible, an estimate of 
the black hole (BH) mass can be obtained by fitting the emitted 


spectrum with a Shakura & S unyaev (1973) accretion disk model 
(Calderone et al. 2013; Sbarrato et al. 2013). 

The large redshift (z > 4) also moves the observed SED 
to lower frequencies. Therefore, the peak of the high-energy 
component appears well below 100 MeV (Ghisellini et al. 
2010a, 2010b; Sbarrato et al. 2012). This introduces a challenge 
in identifying and classifying high-power, high-redshift blazars, 
since a classic hallmark of a blazar is detection by a y- 
ray instrument like the Large Area Telescope (LAT) on the 
Fermi satellite (Atwood et al. 2009). Since the high-energy 
component peaks below 100 MeV, the identification of powerful 
high-redshift blazars in y-rays is problematic. 

On the other hand, sensitive hard X-ray telescopes can detect 
the high-energy hump of high-redshift, extremely powerful 
blazars. In fact, the Burst Alert Telescope (BAT) on board 
the Swift satellite (Gehrels et al. 2004) has detected blazars 
up to larger redshifts than the LAT (Ajello et al. 2009). The 
BAT blazars have hard X-ray spectra (i.e., ax ~ 0.5, assuming 
F(v) oc v~ a A, and this, together with a strong X-ray to optical 
flux ratio, can be generally taken as a signature of the blazar 
nature of a source. 

The most distant blazar known is Q0906+6930 (Romani et al. 
2004; Romani 2006), located at z = 5.47. It was first classified 
as a blazar through a serendipitous EGRET 3a detection. 
However, the more sensitive Fermi /LAT instrument did not 
detect the source, so the EGRET detection could either be 
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spurious or due to an episode of exceptional activity. Subsequent 
to its classification, Q0906+6930 was confirmed as a blazar 
through its X-ray and radio features (Romani 2006). In addition 
to being radio-loud, this source indeed shows a hard X-ray 
spectral index and high X-ray flux, leading to its classification 
as the most distant known blazar. 

We identified the second most distant blazar known, B2 
1023+25, at z = 5.3 (Sbarrato et al. 2012, hereafter S12). 
Its extremely large mass, M B h = 2.8 x 10 9 M 0 , derived by 
fitting the accretion disk spectrum (see SI 2), makes this object 
particularly interesting, since it is possible to put relevant 
constraints on supermassive BH formation models using B2 
1023+25 as a tracer of a population of very high-redshift, 
extremely massive BHs. Indeed, the observation of a single 
blazar with viewing angle 0 V smaller than or comparable with 
the jet beaming angle (i.e., 9 W < 1 /L. where F is the bulk Lorentz 
factor of the relativistic jet) indicates the possible presence of 
2r 2 analogous radio-loud, extremely massive AGNs with their 
jets directed in random directions. A typical blazar has T ~ 15, 
so finding even a few blazars at very high redshift is statistically 
very important for studying the population of extremely massive 
BHs. Therefore, this line of research could become competitive 
with searches for supermassive BHs at high redshifts using 
radio-quiet quasar samples: each blazar with M BH > 10 9 M q 
implies the presence of hundreds of analogous BHs in systems 
with a jet pointing elsewhere. Note that the usual radio-loud 
to radio-quiet ratio (10%) refers to objects with any BH mass. 
At the high-mass end, and at high redshift, this ratio could 
be larger (Volonteri et al. 2011; Ghisellini et al. 2013). The 
existence of z > 5 massive (M BH > 10 9 M G ) BHs in sources 
with powerful jets also raises challenges for BH growth models. 
Rapidly spinning (Kerr) BHs are often invoked as the energy 
source behind powerful jets (e.g., Wilson & Colbert 1995; Sikora 
et al. 2007). However, their accretion disks are radiatively very 
efficient (id = i]Mc 2 , with;; > 0.1, up to 0.3; Thorne 1974). As 
a consequence, over the same time interval, Eddington-limited 
Kerr BHs accrete less matter than Schwarzschild BHs. If this 
is the case, the time needed to form a M BB ^ 10 9 M q BH in 
a source with powerful jets is very long and no such systems 
should exist at z > 4 if they grow primarily through persistent 
accretion (Ghisellini et al. 2013). 

B2 1023+25 was selected from a sample of z > 4 radio- 
loud sources as the best blazar candidate (Sbarrato et al. 2013; 
note that it was also listed as one of the most radio-loud 
quasars at z > 4 by Wu et al. 2013). We were able to observe 
the rising part of its high-energy hump due to a target of 
opportunity (ToO) observation with the X-Ray Telescope (XRT; 
Burrows et al. 2005) on board the Swift satellite (Gehrels et al. 
2004). The hardness of the X-ray spectrum suggested that B2 
1023+25 is a bona fide blazar. However, Swift/X RT observes 
at frequencies too low to properly sample the high-energy 
hump. This motivated our observation of B2 1023+25 with 
the Nuclear Spectroscopic Telescope Array (NuSTAR; Harrison 
et al. 2013). Indeed, a broader X-ray bandpass is required 
to properly classify B2 1023+25 by observing closer to the 
peak of the IC component. Due to its unparalleled broad-band 
sensitivity, NuSTAR is an ideal instrument to determine if the 
X-ray spectrum and flux of high-redshift candidates are hard 
and intense enough, respectively, to classify them as powerful 
blazars. 

Here, we present NuSTAR observations of B2 1023+25 along 
with simultaneous observations obtained in multiple energy 
bands: X-ray observations from Swift/X RT, radio observations 


at three different frequencies from the Combined Array for 
Research in Millimeter-wave Astronomy (CARMA) and the 
Owens Valley Radio Observatory (OVRO), and seven-band 
optical- near-infrared (NIR) photometry from La Silla Observa- 
tory in Chile with the Gamma-Ray Burst Optical Near-Infrared 
Detector (GROND; Section 2). The X-ray data in particular 
allow us to constrain the high-energy component, thereby pro- 
viding important insights into the orientation and Compton- 
boosting of B2 1023+25 (Section 3). 

In this work, we adopt a flat cosmology with Hq = 70 km s -1 
Mpc” 1 and £2 m = 0.3. 

2. OBSERVATIONS AND DATA ANALYSIS 

We performed simultaneous observations of B2 1023+25 with 
five instruments at different frequencies on UT 2013 January 1. 
The X-ray band was covered with Swift/X RT to check for 
possible variability with respect to our previous observations 
and NuSTAR to study the hard X-ray energy range. We added 
data from a previous Chandra observation (Wu et al. 2013) in 
order to increase the statistics of the soft X-ray energy band. 
We re-observed the source from La Silla (Chile) with GROND 
(Greiner et al. 2008), which provides simultaneous photometric 
data from the NIR to the optical in seven different bands. This re- 
observation was performed to check for possible variability of 
the thermal emission of the source. To have simultaneous radio 
data at three different frequencies, we observed B2 1023+25 
with CARMA (Bock et al. 2006) and with the 40 m telescope 
at the OVRO. CARMA observed at 31 and 91 GHz (1 cm 
and 3 mm), while OVRO provided data at 15 GHz (2 cm). 
These data are combined with data from the Wide-field Infrared 
Survey Explorer (WISE 11 ) satellite (Wright et al. 2010) and with 
archival data from NASA/IPAC Extragalactic Database and the 
ASI Science Data Center (ASDC 18 ). 

2.1. X-ray Observations 
2.1.1. NuSTAR Observations 

The NuSTAR satellite observed B2 1023+25 beginning on 
UT 2012 December 31 (sequence 60001107002) for a net 
exposure time of 59.3 ks. The two datasets obtained with 
the NuSTAR Focal Plane Modules A and B (FPMA and 
FPMB) were first processed with the NuSTARDAS soft- 
ware package (v. 1.2.0), jointly developed by the ASDC and 
the California Institute of Technology (Caltech). Event files 
were calibrated and cleaned with standard filtering criteria with 
the nupipeline task using version 20130509 of the NuSTAR 
Calibration Database (CALDB). 

The FPMA and FPMB spectra were extracted from the 
cleaned event files using a circular aperture with a 12 pixel 
(~30") radius, while the background was extracted from two 
distinct nearby circular regions with 30 pixel radii. The ancillary 
response files were generated with the numkarf task, applying 
corrections for the point spread function (PSF) losses, exposure 
maps, and vignetting. The source was detected up to 20 keV 
and the source spectrum in the 4-20 keV energy band was 
formed from a total of 79 counts (of which ~44 are from the 
background) for FPMA and 113 counts (of which ~58 are from 
the background) for FPMB. Both spectra were binned to ensure 
a minimum of 1 count per bin. 


17 Data retrieved from the WISE All-Sky Source Catalog: 
http://irsa.ipac.caltech.edu/. 

18 http://tools.asdc.asi.it/ 


2 


The Astrophysical Journal, 777:147 (8pp), 2013 November 10 


Sbarrato et al. 


2.1.2. Swift Observations 

The Swift satellite observed the source three times: on UT 
2012 June 21 (sequence 00032500001), on UT 2012 June 
22 (sequence 00032500002), and on UT 2012 December 31 
(sequence 00080499001). All XRT observations were carried 
out using the most sensitive photon counting readout mode. 

The XRT dataset was first processed with the XRTDAS soft- 
ware package (v.2.8.0) developed at the ASDC and distributed 
by HEASARC within the HEASoft package (v. 6.13). Event 
files were calibrated and cleaned with standard filtering crite- 
ria with the xrtpipeline task using the latest calibration files 
available in the Swift CALDB. 

The spectra obtained from the single observations are per- 
fectly consistent, with an uncertainty on each measurement of 
~20%-25%, showing no variability among the three observa- 
tions. We therefore merge the individual XRT event files, using 
the XSELECT package for a total net exposure time of 20.3 
ks. Next, we extracted the average spectrum from the summed 
cleaned event file. Events for the spectral analysis were se- 
lected within a circular aperture with a 10 pixel (~23 ff ) radius, 
which encloses about 80% of the PSF, centered on the source 
position. The background was extracted from a nearby circu- 
lar region with a 100 pixel radius. The ancillary response files 
were generated with the xrtmkarf task applying corrections for 
the PSF losses and CCD defects using the cumulative exposure 
map. The latest response matrices available in the Swift CALDB 
were used. The source spectrum in the 0.3-10 keV energy band 
was formed from a total of 41 counts (of which ~3 are from the 
background) and it was binned to ensure a minimum of 1 count 
per bin. 

2.1.3. Chandra Observations 

B2 1023+25 was observed by Chandra on 201 1 March for a 
total of ~5 ks with the ACIS camera. These data were presented 
in Wu et al. (2013). In order to use them together with our other 
data sets, we re-extracted the Chandra spectrum. The data were 
reduced with the CIAO 4.4 package (Fruscione et al. 2006) 
using the Chandra CALDB version 4.4.7, adopting standard 
procedures. The source spectrum was extracted in a circular 
region centered on the peak of the X-ray source emission with a 
radius of 3". The background spectrum was extracted from four 
circular regions with ~5"radii, located around the source. The 
source spectrum in the 0.5-7 keV energy band was formed by 
a total of 54 counts (of which ~ 1 is from the background) and 
it was binned to ensure a minimum of 1 count per bin. 

2.1.4. X-ray Spectral Analysis 

A comparison of the current and previous Swift /XRT ob- 
servations discussed in S12 show that the source did not vary 
between these epochs. If we fit the data of each satellite alone, 
we find a good agreement among NuSTAR, XRT, and Chan- 
dra, but due to the faintness of the source the uncertainties are 
quite large. Since there is no evidence for variability, we per- 
formed a simultaneous fit of the Swift /XRT, Chandra, and NuS- 
TAR spectra using the XSPEC package and adopting C-statistics 
(Cash 1979). 

In prior work, the XRT and Chandra data were fit with a 
simple power-law model plus Galactic absorption (S12; Wu 
et al. 2013). In the present analysis, if we leave A H free to 
vary, we find a high value for the absorption and a slightly 
steeper spectrum than previously reported. The statistical quality 
of the data is, however, not sufficient to distinguish between 


Table 1 

Parameters of the X-Ray Spectral Analysis 


A h 

(cm -2 ) 

T'norm 

(photons cm -2 s -1 ) 

r x 

•^5-10 kev 

(erg cm -2 s -1 ) 

X 2 /dof 

1.5 x 10 20 fixed 
2.8+ 2 j° 7 x 10 21 

L29 -0 2 26 X 10_5 
2.262q77 x 10_5 

1 29 +014 
l - Zy -0.15 

1 60 +0-27 
i -ou_ 0 26 

5.8 x nr 14 
5.5 x to- 14 

230.5/253 

211.3/252 


Note. The errors are at a 90% level of confidence for one parameter of interest. 


no absorption and spectral curvature and a higher level of 
intervening absorption and a steeper spectral index. X-ray 
absorption in the host of blazars is unlikely (the jet is able to 
completely ionize the host interstellar medium), but intervening 
material in high-redshift objects could be responsible for extra 
absorption. The presence or absence of extra absorption due 
to intervening material in quasar X-ray spectra is a matter of 
debate (see e.g., Vignali et al. 2005; Shemmer et al. 2005, 2006; 
Yuan et al. 2006; for a different point of view, see Behar et al. 
2011), so we investigate fits with the column both fixed to the 
Galactic value and left free to vary. Both models provided a 
good description of the observed spectra (see Table 1). If we 
leave Ah free to vary, we find a high value of 2.8 x 10 21 cm 2 
assuming the absorbing material is at z = 0 (for higher redshift 
absorbers, this column increases). A simple power law plus 
Galactic absorption also provides an acceptable fit to the data: 
X 2 = 230.5 for 253 degrees of freedom (dof), to be compared 
with / 2 = 2 1 1 .3 for 252 dof for the case where Ah is left free to 
vary. Although the x 2 clearly improves when Ah is free to vary, 
the reduced / 2 is already < 1 with Ah fixed to the Galactic value; 
therefore, we cannot discriminate between the two possibilities 
from a statistical point of view. The results of the spectral fits 
are shown in Table 1 . The photon spectral index in the two cases 
varies from T x = 1.29^*^ to T x = 1 . 60 *q 26 - following 
analysis, we take this level of uncertainty in the spectral index 
into account, in particular in constraining the jet viewing angle 
(see Section 3). In Figure 1, we plot the X-ray SED of the source 
as derived with the spectral fit performed where Ah is left free 
to vary. 

2.2. GROND Observations 

The 7-band GROND imager, mounted at the 2.2 m MPG/ESO 
telescope at La Silla Observatory, already stated previously 
started observing B2 1023+25 on UT 2013 January 1 at 07:38:43 
UTC. We carried out three 8 minute observations simultaneously 
in all seven g', r’ , f . J, H , andK s bands for a total exposure 
time of 1379 s in the optical and 1440 s in the NIR bands. 
Observations were carried out in an average seeing of 1"2, 
evaluated from the /-'-band image, and at an average airmass 
of 1.8. The source was clearly detected in all bands but g' , for 
which an upper limit of 23.4 (AB magnitude) was found. 

The GROND optical and NIR image reduction and photom- 
etry were performed using standard IRAF tasks (Tody 1993), 
similar to the procedure described in Kriihler et al. (2008). A 
general model for the PSF of each image was constructed using 
bright field stars and it was then fit to the point source. The abso- 
lute calibration of the g' , r' , i' , and z' bands was obtained with 
respect to the magnitudes of Sloan Digital Sky Survey (SDSS) 
stars within the blazar field while the J, H, and K s band cali- 
bration was obtained with respect to magnitudes of Two Micron 
All Sky Survey stars (Skrutskie et al. 2006). 

Table 2 reports the observed AB magnitudes, not corrected 
for the Galactic extinction of E(B — V) = 0.02 from Schlegel 
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Table 2 

GROND AB Observed Magnitudes of B2 1023+25, Obtained UT 2013 January 1 



g' 

r' 

V 

t! 

J 

H 

K s 

^eff (A) 

4587 

6220 

7641 

8999 

12399 

16468 

21706 

AB magnitude 

>23.4 

22.16 + 0.16 

19.91 ± 0.04 

19.73 ± 0.03 

19.52 ± 0.06 

19.20 ± 0.09 

19.35 + 0.24 


Notes. The first row lists the effective wavelengths of each filter (in angstroms). Magnitudes not corrected for Galactic extinction. 



Figure 1. X-ray spectrum of B2 1023+25, along with the two SED models 
discussed in the text. NuSTAR /FPMB data are filled circles (red, circled in 
blue in the electronic version), Chandra data are empty diamonds (blue in the 
electronic version), and Swift/XRT data are open squares. The solid (blue) line 
is the model with 0 V = 3°, V = 13 and parameters as in the first row of Table 3. 
The dashed (red) line is the model with 6 V = 8°, T = 10 and parameters as in 
the second row of Table 3. 

(A color version of this figure is available in the online journal.) 


et al. (1998). Note that these data are fully consistent with those 
from S12, showing that the thermal emission of the source did 
not vary between the two observing times, as was also found for 
the X-ray non-thermal emission (see Section 2.1.4). 

2.3. CARMA Observations 

We observed B2 1023+25 at 31 and 91 GHz (1 cm and 
3 mm, respectively) with CARMA. The observations were 
carried out simultaneously with the NuSTAR observation on 
UT 2013 January 1, with the array in the SL configuration. This 
configuration includes eight 3.5 m antennae on baselines of 5 to 
85 m. Single sideband receivers were used to observe the upper 
and the lower sideband at 3 mm and 1 cm, respectively, and 
the correlator was configured to process 8 GHz of bandwidth. 
After flagging the bad data intervals, the total observation 
time was 3 hr in each band. Strong nearby sources 0956+252, 
0927+390, and the planet Jupiter were used as gain, passband, 
and flux calibrators, respectively. The data were processed using 
the Multichannel Image Reconstruction Image Analysis and 
Display (Sault et al. 1995) software, optimized for CARMA. The 
observations reached an rms of 0.7 (1.5) mJy in the 1 cm (3 mm) 
band, providing a 47a (9 a) detection of the target. The absolute 
flux density calibration, however, adds a systematic uncertainty 
of 10%, so the flux density values used in further analysis are 
/„(31 GHz) = 33 ± 4 mJy and /„(91 GHz) = 14 ± 3 mJy. 


2.4. OVRO 40 m Observations 

The OVRO 40 m telescope obtained a 15 GHz observation of 
B2 1023+25 simultaneous with NuSTAR on UT 2013 January 1. 
The telescope uses off-axis dual-beam optics and a cryogenic 
high electron mobility transistor low-noise amplifier with a 
3 GHz bandwidth. The two sky beams are Dicke switched 
using the off- source beam as a reference and the source is 
alternated between the two beams in an ON-ON fashion to 
remove atmospheric and ground contamination. A noise level of 
approximately 3-4 mJy in quadrature with about 2% additional 
uncertainty mostly due to pointing errors is achieved in a 70 s 
integration cycle. The weighted average of nine consecutive 
integrations was used to derive the 15 GHz flux density 
/„( 1 5 GHz) = 55 ± 4 mJy, where the systematic uncertainty 
in the absolute flux calibration has already been included. 
Calibration is routinely achieved using a temperature- stable 
diode noise source to remove receiver gain drifts and the flux 
density scale was derived from observations of 3C 286 assuming 
the Baars et al. (1977) value of 3.44 Jy at 15 GHz. Details of the 
reduction and calibration procedure can be found in Richards 
et al. (2011). 

3. DISCUSSION 

Figure 1 shows the X-ray data of B2 1023+25. The X-ray 
SED data points were all absorption corrected and rebinned 
to have a 3a detection in each bin. Note that the Swift/XRT 
data are similar to those reported in S12 since variability is 
negligible. Given the rapid variability seen routinely in blazars 
at many wavelengths, this lack of variability is evidence that the 
electrons responsible for the hard X-ray emission have relatively 
small energies and thus lose energy slowly. This is consistent 
with X-rays produced through the so-called External Compton 
process (Sikora et al. 1994), in which relatively cold electrons 
scatter broadline photons. Also, the optical-UV emission is 
steady (see Section 2.2) for a completely different reason. This 
radiation is, in fact, emitted by the accretion disk, which is not 
expected to vary on short timescales. 

The radio part of the spectrum shows flux variability both 
at 5 GHz and at high frequencies (see Figure 2). A 43 GHz 
flux density measurement was published recently by Frey 
et al. (2013), based on a Very Large Array A-configuration 
observation on UT 2002 June 19. A comparison of their 
f y (43 GHz) = 55 ± 4 mJy measurement with our CARMA 
measurements clearly shows that the radio flux is variable in 
time, as expected in blazars. The three radio data points obtained 
in this work define a spectral index a r ~0.7 (F(v) oc v~ ar ), 
steeper than the a r ~ 0.4 value reported by Frey et al. (2013). 
At least in part, this could be due to the fact that the observed 
frequency of 9 1 GHz corresponds to ~570 GHz in the source rest 
frame, likely sampling the optically thin part of the synchrotron 
spectrum. 

Due to the new X-ray flux results, we confirmed the ex- 
treme radio-loudness of B2 1023+25. During the sample 
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Table 3 

Input Parameters Used to Model the SED 


r 

(i) 

6> v 

(2) 

^diss 

(3) 

Rdiss/ RS 

(4) 

p; 

(5) 

B 

(6) 

Vb 

(7) 

Kraax 

(8) 

Si 

(9) 

S2 

(10) 

log A 
(11) 

log Pb 
(12) 

log A 
(13) 

log Tp 
(14) 

13 

3 

504 

600 

0.01 

2.3 

70 

4e3 

0 

2.6 

45.70 

45.93 

44.16 

46.61 

10 

S 

420 

500 

0.23 

4.4 

20 

4e3 

-1 

2.6 

46.72 

46.11 

45.72 

48.26 

5 

20 

588 

700 

7.0 

5.5 

2e3 

4e3 

-1 

2.6 

47.42 

45.98 

45.28 

47.73 


Notes. Column 1: bulk Lorentz factor, Column 2: viewing angle (degrees). Column 3: distance of the blob from the black hole in units of 10 15 cm. 
Column 4: Rdiss in units of the Schwarzschild radius. Column 5: power injected in the blob calculated in the comoving frame, in units of 10 45 erg s -1 . 
Column 6: magnetic field in Gauss, Columns 7 and 8: minimum and maximum random Lorentz factors of the injected electrons, Columns 9 and 10: 
slopes of the injected electron distribution [<2(yXI below and above yb, Column 11: logarithm of the jet power in the form of radiation, Column 12: 
Poynting flux, Column 13: bulk motion of electrons, and Column 14: protons (assuming one cold proton per emitting electron). The spectral shape of 
the corona is assumed to be oc v -1 exp(— /zv/150keV). For all models, we have assumed a radius of the broadline region of Rblr = 9.2 x 10 17 cm, a 
BH mass of 2.8 x 10 9 M©, and an accretion disk luminosity of Ld = 9 x 10 46 erg s -1 , corresponding to Ld/L^dd = 0.25. 



Figure 2. Broad-band SED of B2 1023+25 together with the models discussed 
in the text. Simultaneous OVRO, CARMA, GROND, and NuSTAR data are 
filled circles (red points circled in blue in the electronic version). Chandra data 
are open diamonds (blue in the electronic version), while Swift/XKT data are 
open squares. The (gray) filled symbols are data from the literature: squares 
are archival data from ASDC, the diamond is the radio point from Frey et al. 
(2013), the circles and the two upper limits are WISE data, and the line is the 
SDSS spectrum. The dotted (black) line is the thermal emission of the source, 
including the accretion disk, torus, and X-ray corona emission. The Fermi /LAT 
upper limit is for 3.8 yr, 5 o (red arrow). The solid (blue) line is the model with 
parameters as in the first row of Table 3. The dashed (red) line is the model with 
parameters as in the second row of Table 3. The dot-dashed (green) line is the 
model with parameters as in the third row of Table 3. 

(A color version of this figure is available in the online journal.) 


selection (Sbarrato et al. 2013), we used the canonical 
radio-to-optical ratio to define its radio-loudness ( R = 
A [5 GHz]/ A [2500 A] ~ 5200) and this allowed to classify 
B2 1023+25 as the most radio-loud quasar of our sample. 
In addition to this, we now calculate the X-ray based radio- 
loudness R x = vL v [5 GHz]/ Lx [2-10 keV], using the X-ray 
fluxes and spectral indices listed in Table 1 . We obtain log R x = 
—0.65 and log R x = —0.72 (for fixed and free Ah, respec- 
tively). Both values confirm the extreme radio-loudness of B2 
1023+25 according to the calibration introduced by Terashima & 
Wilson (2003), which classifies quasars as radio-loud if they 
have log R x > —4.5. 

In S12, we derived a set of parameters that reproduced the 
observed SED and suggested classifying B2 1023+25 as a blazar 
(bulk Lorentz factor F — 14; jet viewing angle (9 V = 3°). We 


fit the new observations using the model described in Ghisellini 
& Tavecchio (2009). Since this is a one-zone model, which 
assumes that the emitting region is rather compact, it cannot 
account for radio emission, which in the considered region is 
severely self-absorbed. 

In this model, both Q y and T are free parameters and can be 
chosen independently. 

Because of the hard and bright X-ray spectrum shown by 
B2 1023+25, we find a small value of 9 V and a large Doppler 
boosting (i.e., large T). We find 6 V < 1/T, as is typical of 
known blazars. Hence, we confirm B2 1023+25 as a blazar (see 
Section 3.1). 

Because of the limited statistics in the X-ray spectrum, we 
investigate the range of models consistent with the uncertainties. 
As noted above, depending on how the spectrum is modeled, 
the intrinsic continuum may be softer and fainter overall. This 
case implies a larger value of 6 V and a somewhat smaller value 
of r (see Section 3.2). The jet viewing angle, 0 V , associated with 
this limiting solution is an upper limit. Since this model is also 
characterized by less Doppler boosting, it corresponds to a larger 
intrinsic luminosity relative to the SED corresponding to the 
X-ray best fit parameters. We consider then “re-orienting” the 
jet to a typical blazar viewing angle (i.e., ~3°) and we check if 
the corresponding SED resembles the one of a typical powerful 
blazar seen at lower redshift. We then use this to check the 
reliability of the obtained solution; that is, we require that, if the 
jet was pointed toward us, the solution would show reasonably 
similar properties to the blazar sample. 

In our modeling, we keep the parameters associated with 
the thermal emission from the accretion disk fixed. We assume 
a BH mass Mbh = 2.8 x 10 9 M Q and an accretion disk 
luminosity L d = 9 x 10 46 ergs -1 , as derived in S12. Note that 
varying the BH mass value inside the formal confidence range 
(Mbh = 1.8-4. 5 x 10 9 M 0 ) does not change the results of our 
SED modeling. 

3.1. Best Fit: Small Viewing Angle, Large Bulk Lorentz Factor 

In our best fit model, we find a set of parameters consistent 
with those from S12 (T = 13, 6 V = 3°). We report these in the 
first line of Table 3 as the best fit to the broad-band SED. The 
case in Table 3 corresponds to the best fit to the X-ray data with 
Ah left free to vary. Using Ah fixed to the Galactic value yields 
a harder spectrum and therefore an even more extreme blazar 
classification. The model (blue solid line in Figures 1, 2, and 3) 
describes a typical blazar, with a viewing angle smaller than the 
jet beaming angle (0 V < 1/T), firmly classifying B2 1023+25 
as a blazar. 
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Figure 3. SED of B2 1023+25 with data as in Figure 2. The curved (gray) stripe 
corresponds to the sensitivity of Fermi /LAT after 1 yr of operation (5a, lower 
bound) and 3 months (10a, upper bounds). The solid (blue) line is the model 
with parameters as in the first row of Table 3. The dashed (red) line is the model 
with parameters as in the second row of Table 3 (i.e., 6 V = 8°), but “re-oriented” 
at 6 W = 3°, as labeled. The dot-dashed (green) line is the model with parameters 
as in the third line of Table 3 (i.e., 0 V = 20°), but “re-oriented” at 0 V = 3°, as 
labeled. 

(A color version of this figure is available in the online journal.) 

If 9 V < 1 /r, as in our best-fit model, the number of sources 
similar to B2 1023+25 but with the jet oriented outside of our line 
of sight is 2r 2 = 338(F / 13) 2 in the portion of the sky covered 
by SDSS+FIRST (Ghisellini et al. 2010a; Volonteri et al. 2011; 
Ghisellini et al. 2013). Since the combined SDSS and FIRST 
surveys (from which we selected B2 1023+25) cover together 
8770 deg 2 , this implies that in the whole sky there must be at 
least ~ 1550 sources that share the same intrinsic properties of 
B2 1023+25. Since the co-moving volume in the redshift frame 
5 < z < 6 is ~ 3 80 Gpc 3 , we can conclude that there must 
be at least four radio-loud AGNs similar to B2 1023+25 per 
Gpc 3 . Albeit extrapolating from a sample of one, this would 
imply the presence in the redshift bin 5 < z < 6 of at least four 
supermassive BHs per Gpc 3 , with a BH mass of M BH ~ 10 9 M 0 , 
hosted in jetted systems. 

3.2. Slightly Misaligned Jet 

Figure 1 shows that the X-ray data have large error bars. A 
softer X-ray spectrum cannot be excluded at a 90% confidence 
level (see Table 1). A softer spectrum implies a larger viewing 
angle and therefore a somewhat less extreme bulk Lorentz 
factor. Specifically, a viewing angle of 6 V = 8° with T — 
10, together with the other parameters in the second line 
of Table 3, still reproduce the broad-band data. The model 
(red dashed line in Figures 1 and 2) represents an alternative 
interpretation consistent with the X-ray data points at a 90% 
level of confidence. This viewing angle, slightly larger than 
the jet beaming angle 1/T, still classifies B2 1023+25 as a 
blazar. As a consistency check, we test how an object with the 
same intrinsic (comoving) properties would look if oriented 
at 0y = 3°, i.e., at 0 V < 1/T. The re-oriented model is 
shown in Figure 3 by the dashed line (red dashed line in the 
electronic version). The resulting X-ray flux would be unusual 
but not unprecedented, being very similar, for example, to GB 
1428+4217 ( z = 4.72; Worsley et al. 2004), although the latter 
shows a synchrotron hump much dimmer than our “re-oriented” 
B2 1023+25. We conclude that 0 V = 8° is the largest possible 
viewing angle consistent with the NuSTAR data. 


If 0y > 1/T, as in the above case, the number of sources 
similar to B2 1023+25 but with their jet oriented in random 
directions is 1/(1 — cos Of). Hence, with 9 V = 8°, in the portion 
of sky covered by SDSS+FIRST, there would be 103 AGNs 
analogous to B2 1023+25 and in the whole sky there would 
be 469 (i.e., ~1.2 object per Gpc 3 ). Even in this limiting case, 
the number of extremely massive BHs in jetted systems in the 
redshift bin 5 < z < 6 is cosmologically significant. 

3.3. Can the Jet be at 20° from our Line of Sight? 

Frey et al. (2013) claim that a viewing angle of at least 
~20° with a bulk Lorentz factor F ~ 15 can be inferred 
for B2 1023+25 from published 5 GHz very long baseline 
interferometry (VLBI) imaging data. We therefore attempt to 
fit our X-ray data with 6 V = 20° to test this hypothesis. We find 
that the data are not consistent with both a large viewing angle 
and a large Lorentz factor, i.e., the values obtained by Frey et al. 
(2013). 

In this case, the corresponding Doppler factor S = [F( 1 — 
ft cos 0y )\~ 1 ~ 1 and the intrinsic jet power become huge, to 
account for the observed X-ray flux (i.e. , Pjet ~ 10 50 erg s *). 
In addition, the fit to the observed data is poor. 

The maximum bulk Lorentz factor providing a good fit to the 
NuSTAR and broad-band data with a viewing angle of 0 V = 20° 
is T = 5 (along with the parameters in the third line of Table 3). 
This model is shown in Figure 2 as the dot-dashed line (green 
dot-dashed line in the electronic version). The corresponding 
beaming factor is 5 ~ 2.5. Such a modest beaming factor implies 
that the intrinsic luminosity would be very high. This would 
imply a class of objects with an extreme intrinsic luminosity. 
If such objects existed (at any redshift), we should see a few 
of them pointing at us. For illustration, we then “re-orient” B2 
1023+25 to 0y = 3° (see Figure 3). Similar SEDs have never 
been observed, at any redshift. All powerful blazars observed so 
far have the Compton component dominating the overall SED, 
contrary to what is shown in Figure 3. We therefore believe that 
it is highly unlikely that B2 1023+25 can be described with 
0 V = 20° in the high-energy emitting region. We cannot exclude 
the possibility that the jet bends between the X-ray and the radio 
emitting regions. In this case, it is possible that the large-scale 
jet (i.e., the radio emission) is seen at a larger viewing angle 
than the compact jet. 

Furthermore, consider that the 5 GHz VLBI observations 
analyzed by Frey et al. (2013) correspond to a rest-frame 
frequency of 31.5 GHz. At this frequency, all of the VLBI 
components except for the very inner core are emitting thin 
synchrotron radiation. Since the brightness temperature of a 
synchrotron source peaks at the self-absorption frequency, we 
conclude that all the brightness temperatures of the resolved 
components are lower limits. For the core, Frey et al. (2013) 
indeed performed a fit with a resolved plus an unresolved 
component. It is very likely that the resolved core is optically 
thin (thus giving a lower limit to the brightness temperature), 
while the unresolved core gives a lower limit because of the 
upper limit on the size. As a consequence, the derived Doppler 
factors are all lower limits and the derived viewing angles are 
all upper limits. 

4. CONCLUSIONS 

We selected B2 1023+25 as the best z > 4 blazar candidate 
from the SDSS+FIRST quasar catalog and we classified it as a 
blazar as the result of a Swift /XRT ToO observation (S12). 
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Here, we report simultaneous NuSTAR and Swift observations 
to improve the broad-band X-ray spectrum and further cement 
the blazar classification. We use the improved data to determine 
the jet orientation and the relativistic boosting factor. Simultane- 
ous GROND data are important to check for possible variability. 
From the comparison with our first data, B2 1023+25 does not 
show variability in either its thermal or non-thermal emission. 

We fit the broad-band SED with the model described in 
Ghisellini & Tavecchio (2009), focusing on the X-ray energy 
band as an important constraint, and we analyze in detail the 
non-thermal jet emission of the source. We confirm that B2 
1023+25 is an extremely radio-loud quasar, with a jet oriented 
very close to our line of sight, and hence the Doppler boost 
is large. Our SED modeling indicates a small viewing angle 
(f) v = 3°) associated with a large bulk Lorentz factor of r = 13. 
To account for the large X-ray data uncertainties, we tested 
solutions with larger viewing angles. A model with 0 V = 8° and 
T = 10 cannot be excluded at the 90% uncertainty level by the 
data. A viewing angle larger than this is not consistent with the 
data and the resulting solution provides a lower limit to the real 
X-ray spectrum. Therefore, B2 1023+25 shows a jet orientation 
and a Doppler boosting that lead us to firmly classify it as 
the second most distant blazar known (z = 5.3). This implies 
the presence in the SDSS+ Faint Images of the Radio Sky at 
Twenty-centimeters (FIRST) survey of several detectable radio- 
loud sources with jets oriented in other directions. However, we 
have not been able to identify such objects in the survey; there 
are only four other radio-detected quasars in the SDSS+FIRST 
sample at z > 5. Although the statistics are small and possibly 
not constraining, the apparent inconsistency is addressed in 
Volonteri et al. (201 1). 

B2 1023+25 is the first object at z > 5 detected by NuSTAR 
and confirms that NuSTAR is a very useful instrument to deepen 
our knowledge of the high-redshift X-ray universe. Specifically, 
NuSTAR could be an ideal tool to continue the z > 4 blazar 
hunt. With the study of this single object, indeed, we were able 
to estimate, albeit with large uncertainty, how many extremely 
massive BHs in jetted sources are present at 5 < z < 6. This 
constrains the mass function of heavy BHs injetted systems as a 
function of redshift, which provides a complementary constraint 
for surveys of radio-quiet AGNs (see e.g., Ghisellini et al. 2013; 
Willott et al. 2010). The confirmation that B2 1023+25 is a 
blazar strengthens the suggestion of Ghisellini et al. (2013) that 
there are two epochs of heavy BH formation: radio-loud objects 
preferentially form their M BH > 10 9 M o BHs at z ~ 4, while 
radio-quiet quasars are formed at z ~ 2. 
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